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Abstract 

(N 

^ Jn this paper we discuss some simple analytical models to fit the cosmic-ray (CR) proton data collected 
D by the AMS detector in June 1998 and by the PAMELA detector in several campaigns covering 
the period 2006-2009. The CR proton spectrum at Earth is derived starting from the model of the 
■local interstellar spectrum (LIS) and folding it with the solar modulation potential in the force field 
approximation. The data are well described by a LIS modeled with a simple power law particle 
p~j^ .momentum density. 
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1. Introduction 

Cosmic rays (CRs) interact with gas atoms 
during their propagation in the interstellar 
medium and can suffer significant energy losses, 
!thus modifying their injection spectra and 
composition. In addition, the spectra of CRs 
teaching the Earth are affected by the solar 
wind and the by the solar magnetic field (solar 
modulation effect). The solar modulation plays 
!a relevant role on CR spectra in the low energy 
^region, and its effect needs to be disentangled 
to allow a comprehensive picture to emerge. In 
.fact, to understand the origin and propagation of 
CRs, a knowledge of their energy spectra in the 
interstellar medium is required. 

Precise measurements of CR spectra over a 
wide rigidity range, from a few hundred MV to 
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tens of GV can be used to study the effect of 
solar modulation, including the convective and 
adiabatic cooling effect of the expanding solar 
wind and the diffusive and particle drift effects of 
the turbulent heliospheric magnetic field (HMF). 

A full three-dimensional (3D) model was 
developed to compute the differential intensity of 
CR protons from 10 MeV to 30 GeV at Earth [l|, 
and was applied to give an interpretation of the 
PAMELA proton data sets collected from 2006 
to 2009 jl]. The model also includes a detailed 
treatment of the CR propagation in the solar 
magnetic field, that allows a precise description of 
the solar modulation effect. The implementation 
of this approach requires to provide the local 
interstellar proton spectrum (LIS) as initial 
condition. The LIS input spectrum is then 
"modulated" by the solar magnetic field, that 
affects the shape of the CR spectrum at Earth. 

The choice of the LIS has always been 
rather contentious (see for instance [1]), and its 
parametrization in terms of proton kinetic energy 
could be complex (see for instance [l|, 0], jEf). 
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In this paper we assume some simple analytical 
LIS models to fit the proton fluxes measured 
by the AMS detector in June 1998 0] and by 
the PAMELA detector in different periods from 
2006 to 2009. The solar modulation effect is also 
described in a simple form, using the force-field 
approximation 01 . 

2. Local proton spectrum models 

The simplest model describing CR 
acceleration is the first-order Fermi mechanism, 
where particles gain energy by diffusing back 
and forth across a shock front while convecting 
downstream. In this framework, the particle 
differential density per unit momentum n{p) is oc 
p~°'. After injection into the interstellar medium 
with spectral index a ^ 2.1 ^ 2.2, characteristic 
of supernova remnant (SNR) shocks, CRs are 
transported in the astronomical environments 
with rigidity-dependent escape lengths, that 
soften their spectra by 6 ^ 0.5 -^0. 6, leaving 
a steady-state CR particle momentum density 
n{p) oc where a = a + 5 ~ 2.7 4- 2.8 j^. 

Therefore, in the present work we will consider 
a simple LIS model with the differential particle 
momentum density in the form: 



n{p) 



The spectral differential intensity in momentum is 
obtained by multiplying n{p) for the factor (Sc/Att, 
where (3c is the particle velocity: 



Jip) = T^n{p) = kl3 (— 
4vr \po 



(2) 



where k = koc/A-K. Hereafter we will assume c = 1 
and we will express both energies and momenta in 
units of GeV. In writing the previous equations 
we introduced a momentum scale po = 1 GeV. 
Therefore ko will be expressed in the same units as 
n{p), i.e. in GeV~^m~^ and k will be expressed in 
the same units as J{p), i.e. in GeV~^m~^s~^sr~^. 

The differential intensity in momentum can be 
converted into a differential intensity in kinetic 
energy taking into account that: 
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Computed LIS Intensity vs time 
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Figure 1: Fit of the PAMELA proton data collected in 
November 2006 (red), December 2007 (green), December 
2008 (blue) and December 2009 (cyan) with a simple 
momentum power law LIS folded with solar modulation 
(the fit was performed up to 175 GeV). The top panel 
shows the fit results superimposed to the the data points. 
The dashed lines show the fitted spectra (Eq. [7]) and the 
continuous lines show the corresponding LIS (Eq. [5]). The 
bottom panel shows the fit residuals. 
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^T{T + 2m) 

where m is the particle rest mass. The CR 
spectrum in the interstellar space as a function 
of the kinetic energy is therefore given by: 



Jlis{T) = k 



T{T + 2m) 



(5) 



The spectral index in kinetic energy, a^, is 
defined as: 



ax = —- 



d\ogJus{T) T + m 



dlogT 



a 



T + 2m 



(6) 
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Figure 2: Fit of the AMS proton data collected in June 
1998 (red) with a simple momentum power law LIS folded 
with solar modulation. The top panel shows the fit results 
superimposed to the the data points. The dashed lines 
show the fitted spectra (Eq. [7]) and the continuous lines 
show the corresponding LIS (Eq. [S]). The bottom panel 
shows the fit residuals vs energy. 



The previous result shows that the CR spectrum 
exhibits a change of curvature with increasing 
kinetic energy [ax — >■ a/2 for T — )■ while 
— ^ « for T — )■ oo). 

The LIS model of Eq. O can be folded with 
the solar modulation effect in the force field 
approximation introducing the solar modulation 
potential field $, and results into a CR spectrum 
at Earth given by: 



JEarth{T) = JlIs{T+^Z/A) 



T{T + 2m) 



Figure 3: Fit of the PAMELA proton data collected in 
November 2006 (red), December 2007 (green), December 
2008 (blue) and December 2009 (cyan) with the LIS of 
Eq. [S] folded with solar modulation (the fit was performed 
up to 175 GeV). The top panel shows the fit results 
superimposed to the the data points. The dashed lines 
show the fitted spectra and the continuous lines show the 
corresponding LIS (Eq. [8]) . The bottom panel shows the 
fit residuals. 



of Ref. [21]. The data at higher energies, from 
48 GeV to 1 TeV have been taken from Ref. j^. 
The fit was performed up to 175 GeV using 
the MINUIT package implemented in the ROOT 



(T + m + ^Z/Ay 



(7) 

where Z and A are the atomic and mass number 
respectively of the given CR species. 

Figure [1] shows the results obtained fitting the 
PAMELA proton data with Eq. [3 The data 
up to 48 GeV have been taken from Table 1 



framework [lOj]. The residuals exhibit very small 
fluctuations, within a few %. It is worth to point 
out that the 2009 data show some large point-to- 
point fluctations, probably due to the fact that 
these data were collected at the end of the 11*'* 
' Solar cycle. Figure |2] shows the results obtained 
fitting the AMS proton data with Eq. [3 The data 
have been taken from Table 3 of Ref. j^. The 
error bars shown in the plots are evaluated by 
adding in quadrature statistical and systematic 
uncertainties. Also in this case the residuals 
exhibit small fluctuations. 
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The results of all the fits are summarized in 
Table [TJ It is interesting to point out that both 
the AMS and PAMELA data sets are well fitted 
by the momentum power law LIS. Moreover, 
the values of the LIS parameters (prefactor and 
spectral index) obtained from the fit of the AMS 
data are consistent with those obtained from the 
fits of the PAMELA data sets taken in 2006, 2007 
and 2008, as one would expect since the LIS is 
time independent. On the other hand, the fit 
of the data collected by PAMELA in 2009 yields 
values that differ significantly from those obtained 
from the other fits. Because of the consistency 
between the AMS and the first three PAMELA 
fits, we decided to combine these results into a 
unique LIS, with a spectral index a = 2.853 and 
a pre-factor of 2.4863 x 10^ GeV-^m-^g-isr-^ 
evaluated from the weighted average of the 
individual fit results. Using these parameters 
for the LIS (Eq. [5]) and leaving only the solar 
modulation potential free, the fits of the AMS and 
of the PAMELA data do not change significantly 
with respect to the ones shown in Figures [T] and 
[21 In particular, the fitted values of the solar 
modulation potential do not exhibit significant 
variations (see Table [1]). 

We have also fitted the data samples using the 
LIS model given in Ref. i.e.: 

{ 707 k e4-64-0.08(logT)2-2.91v^ T <To 
685 k e3-22-2.781ogT-1.5/T > 

(8) 

where Tq = 1.4 GeV. The numerical 
coefficients in the previous equation include units 
of measurement: Jlis{T) is given in units of 
GeV~^m~2s~^sr~^ if T is expressed in GeV. The 
coefficient k is a scale factor {k = 1 wilLreproduce 
exactly the formula in Eq. 13 of Ref. jl|). 

Figure |3] shows the results of the fits performed 
assuming the LIS model in Eq. |8] folded with the 
solar modulation in the force field approximation. 
In this case the residuals show larger fluctuations 
with respect to the simple power law fits and 
the values are higher. The fit results are 
summarized in Table [1] It is also worth to point 
out that, when this LIS model is assumed, the 



solar modulation potential values are smaller with 
respect to those obtained assuming the simple 
momentum power law model. This feature is 
due to the shape of the LIS given in Eq. IH that 
predicts a curvature at low energies. Another 
interesting feature of these fits is the value of the 
prefactor, that in all cases is consistent with 1, 
as expected. We have also performed the fit of 
the AMS and PAMELA data with the LIS model 
of Eq. |8] with a fixed prefactor k = 1. This 
constraint does not worsen the fit, and the the 
solar modulation potential values do not change 
significantly. 

To reproduce the B/C ratio the injection 
spectrum of CR protons is often approximated 
with a power law in the rigidity space with a 
break at a few GV(see for instance jij). A 
possible description of this feature can be given by 
choosing for the LIS momentum density a broken 
power law function (with a discontinuity in the 
first derivative at the break): 
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that result into a differential intensity given by: 
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(10) 



with p = \JT{T + 2m) and k = koc/Air. The 
break momentum pb corresponds to a break 
kinetic energy Tb = ^/pl + rn"^ — m. 

The results of the fits with the LIS of Eq. [10] 
are also shown in Table [H We note that the fits 
yield no evidence of a break, since ai ~ 0^2 within 
the errors for all the data sets. 

3. Conclusions 

The simple power law model of the proton LIS 
folded with the solar modulation in the force field 
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approximation provides a good fit of tlie botli 
AMS and PAMELA proton data. Tliis result 
needs to be investigated witli furtlier analyses, 
for instance by using data from helium and heavy 
nuclei on short time periods, that are not publicly 
available. 

In general, it is worth to point out that a CR 
measurement at Earth does not allow to easily 
reconstruct the LIS, since the Solar modulation 
effect cannot be easily disentangled. However, a 
constraint to the LIS spectrum could be provided 
by a fit of the gamma-ray emissivity of the local 



neutral gas measured by the Fermi LAT |8|, 111 



12 



Acknowledgements 

We are grateful to Charles D. Dermer and 
Andrew W. Strong for the fruitful discussion 
during the preparation of the manuscript and for 
their valuable contribution. 



References 

[1] M. S. Potgieter, E.E. Vos, M. Boezio, N. De 
Simone, V. Di Felice and V. Formato Modulation 
of galactic protons in the heliosphere during 
the unusual solar minimum of 2006 to 2009, 
'http : //arxiv . org/abs/1302 . 1284 

[2] O. Adriani et al., Time dependence of the proton flux 
measured 

by PAMELA during the July 2006 - December 2009 

solar minimum, http: //arxiv. org/abs/1301.4108j 
[3] I.V. Moskalenko,"A;W. Strong, J.F. Ormes and M.S. 

Potgieter, Astrphysical Journal 565 (2002) 280 
[4] W.R. Webber and P.R. Higbie, Journal of 

Geophysical Research 115 (2010) A05102 
[5] I.V. Moskalenko and T.A. Porter, Astrphysical 

Journal 670 (2007) 1472 
[6] J. Alcarez et al. (AMS Collaboration), Physics 

Letters B 490 (2000) 27 
[7] L. J. Gleeson and W. I. Axford, ApJ 154 (1968) 1011 
[8] C. D. Dernier, Physical Review Letters 109 (2012) 

091101 

[9] O. Adriani et al.. Science 332 (2011) 69 



[10] http: / /www. root. cern.ch] 



[11] A. A. Abdo et al, Astrophys. J. 703, 1249 (2009) 
[12] Jean-Marc Casandjian, AIP Conf. Proc. 37 (2012) 
1505, |http : //dx . doi . org/10 . 1063/1 . 4772218] 



5 



Simple power law LIS (Eq.[5} 




AMS Jun 1998 


PAMELA Nov 2006 


PAMELA Dec 2007 


PAMELA Dec 2008 


PAMELA Dec 2009 




2.342 ± 0.193 


2.548 ± 0.067 


2.520 ± 0.065 


2.414 ± 0.064 


2.072 ± 0.097 


a 


2.832 ± 0.024 


2.858 ± 0.007 


2.857 ± 0.007 


2.847 ± 0.007 


2.818 ± 0.011 


«(GoV) 


0.629 ± 0.030 


0.685 ± 0.009 


0.602 ± 0.008 


0.557 ± 0.008 


0.454 ± 0.012 


X^/d.o.f. 


1.4/25 


16.2/87 


15.1/87 


11.1/87 


28.2/87 



Simple power law LIS (Eq.[5) with k = 2.4863 X 10* GeV" ^ m "^s" ^ sr " ^ and a = 2.853 




AMS ,Iun 1998 


PAMELA Nov 2006 


PAMELA Dee 2007 


PAMELA Dec 2008 


PAMELA Dee 2009 


«(GcV) 


0.646 ± 0.010 


0.676 ± 0.003 


0.598 ± 0.003 


0.566 ± 0.003 


0.501 ± 0.005 




2.4/27 


17.8/89 


15.5/89 


12.5/89 


59.3/89 



LIS from ref. |1| (Eq.|8) 




AMS Jun 1998 


PAMELA Nov 2006 


PAMELA Dee 2007 


PAMELA Dec 2008 


PAMELA Dee 2009 


k 


1.062 ± 0.025 


1.044 ± 0.078 


1.029 ± 0.080 


1.015 ± 0.080 


0.995 ± 0.089 


«(GeV) 


0.583 ± 0.018 


0.594 ± 0.005 


0.500 ± 0.005 


0.459 ± 0.005 


0.377 ± 0.007 


X^/d.o.f. 


11.6/26 


56.9/88 


82.3/88 


104.4/88 


54.1/88 








LIS from ref. | IJ (Eq.|S} with k = 


1 






AMS ,Iun 1998 


PAMELA Nov 2006 


PAMELA Dee 2007 


PAMELA Dec 2008 


PAMELA Dee 2009 


*( GeV) 


0.550 ± 0.011 


0.572 ± 0.004 


0.488 ± 0.003 


0.453 ± 0.003 


0.379 ± 0.007 


X^/d.o.f. 


17.9/27 


88.3/89 


95.3/89 


107.6/89 


54.4/89 



Broken power law LIS (Eg . llOt 




AMS Jun 1998 


PAMELA Nov 2000 


PAMELA Dee 2007 


PAMELA Dec 2008 


PAMELA Dec 2009 


fc( GeV^-'m-^s--'sr--') 


134.38 ± 25.18 


157.24 ± 16.30 


155.69 ± 10.14 


153.36 ± 10.42 


132.67 ± 29.38 


ai 


2.705 ± 0.155 


2.942 ± 0.089 


2.888 ± 0.053 


2.844 ± 0.051 


2.908 ± 0.112 


a2 


2.824 ± 0.024 


2.863 ± 0.012 


2.858 ± 0.007 


2.847 ± 0.007 


2.817 ± 0.012 


«(GeV) 


0.566 ± 0.080 


0.722 ± 0.039 


0.615 ± 0.023 


0.556 ± 0.021 


0.484 ± 0.043 


Pb(GeV) 


6.13 ± 0.41 


5.96 ± 2.11 


5.94 ± 0.14 


5.91 ± 0.14 


6.00 ± 0.47 


X^/dof 


1.0/23 


14.9/85 


14.9/85 


11.1/85 


27.6/85 



Table 1: Results of the fits of the AMS and PAMELA proton data with the various LIS models discussed in the text. 
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